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ABSTRACT 

We give a simple derivation of the Virasoro constraints in the Kontsevich model, 
first derived by Witten. We generalize the method to a model of unitary matrices, 
for which we find a new set of Virasoro constraints. Finally we discuss the solution 
for symmetric matrices in an external field. 
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1. Introduction. 



There has recently been considerable progress in the study of matrix models, 
following the remarkable proof by Kontsevich that the intersection numbers (correla- 
tion functions) of two-dimensional topological gravity are generated by a new type of 
matrix model [1]. A short time later, Witten showed [2] that the partition function 
of this model obeys the Virasoro constraints of the one-matrix model [3, 4], thus 
completing the chain of arguments in a proof of the old conjecture that topological 
gravity and matrix models are equivalent [5]. 

However, Witten's lengthy proof depends on a cumbersome diagrammatic expan- 
sion. We present in section 2 a much simpler derivation, which arose out of recent 
work into matrix integrals involving an external field [6] . This approach allows us to 
consider other models; and in section 3 we tackle the case of a unitary matrix in an 
external field [8, 6], deriving a set of Virasoro constraints which describe the contin- 
uum limit of the multicritical unitary matrix models [9, 10]. These results are new, 
but after this work was completed, we learned that another group [11] has found the 
same constraints from the complementary viewpoint of the mKdV flows [10, 12]. 

In section 4 we present the results of some work on symmetric matrices in an 
external field. As far as we know this model has never before been investigated, so we 
stray somewhat from the main subject of this paper and give some general details. 
Unfortunately this solution is incomplete, and in particular we fail to produce for 
symmetric matrices any results generalizing the Virasoro constraints of sections 2 
and 3. 

In the final section we discuss our results, and speculate what these might teach 
us about unitary matrix integrals. 



f See also ref. 7. 

§ We thank E. Witten for informing us of this. 
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2. Virasoro constraints from Schwinger-Dyson equations. 



Consider the integral over N x N hermitian matrices 

Z = JvMexpTr (xM - \gM^ , (2.1) 

which is a function of the N eigenvalues x a of the hermitian matrix In refs. 6, 13, 
it is shown how to evaluate Z using the Schwinger-Dyson equations of motion, 

/ v£l M exp Tr (™ " = ( xT " 9 Wc*) i = ' (2,2) 

which can be recast as a set of N differential equations in terms of the eigenvalues 
d 2 Z | 1 (dZ dz\ _ 1^ ^ 

Full details of the derivation are given in refs. 8, 6. In ref. 13 these equations were 
used to find Z in the spherical approximation, and recently we showed how to extend 
this solution to all orders in N [6].* In this paper we start from these same equations, 
but the analysis will be very different. 

Closely related to Z is the integral 

Z = JvMexpTr {-\AM 2 - ^gM 3 ) , (2.4) 

for a positive definite hermitian matrix A. In his ground-breaking paper [1], Kontse- 
vich showed how the expansion of Z in terms of Feynman diagrams can be interpreted 



If Note that for the purposes of this paper, we define the integral without an iV in the exponent. 
* The reader is advised that an early preprint version of ref. 6 contained some mistakes in the 
introduction. 
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as a cell decomposition of the moduli space of Riemann surfaces [14, 15, 16]. By a 
shift of integration variables, it is easy to see that Z and Z are related by 

A 3 

Z = expTr^-2 xZ, (2.5) 

provided we choose X = Our goal now is to use (2.3) and (2.5) to show that Z 
obeys the Virasoro constraints of the one-matrix model [3, 4]. 

Following Kontsevich, we factor Z in the form 

Z = H( f i a + f i b )- 1/2 Y , (2.6) 

a.b 



where {/x a } are the eigenvalues of A. Then, after a change of variables to A a = /i a — 
Agx a , eqns. (2.3) and (2.5) imply that Y satisfies the set of differential equations 

d^Y + 9Y(^_ZA +T ^__ (BY _dY_ ) +Y (J_ + A)= Q 

d\l + d\ a W to) t£rf-ri\d\ b d\ a ) + Ati) ' 

(2.7) 

where we have defined 

^ ^ + fib 2fc + 1 H 

Next we make good use of Kontsevich's important observation that the perturbative 
expansion of Y depends on the eigenvalues of A only through the invariants t k . This 
lets us change variables from A a to t k , i.e., 



and after some algebra we obtain 

i „ i A i 



2 S , 2fe+3 f)f, ' ( 2 - 9 ) 



4 2fc+2j+6 d k9jY + ^ i 2k+2 9kY 

k,j=0 °y k=0 ^a 

fc=0 x ^ a ^ a 



(2.10) 
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where d k = d/dt^. Note that these equations are formally valid independent of the 
large-iV limit, except that only for N — > oo are the t^s truly independent, and for 
finite N the sums over k must be truncated.** 

To make contact with Kontsevich's results, we set g = i/2. Then we can pick off 
the coefficient of 1/ /i„ +2n to obtain the equations 

oo 

n = -1 : Y^(k + \)t k d k ^Y + \t 2 Q Y = \d Q Y , 
k=i 

oo 

n = o : ^(fc + i^y + ^y^i^y , (2.11) 

00 n 

n > 1 : J2( k + 5)*fc 9 *+« y + 3 S = 5^+1^ ■ 

fe=0 fc=l 

Given that we have used a different normalization for the t^'s, these are precisely the 
equations first derived for the Kontsevich model in ref. 2, 

L n Y = \d n+1 Y, n>-l, (2.12) 

i.e., the Virasoro constraints of refs. 3, 4, corresponding to perturbations about the 
m = 1 "topological" point of the one-matrix model. It is remarkable that the entire 
set of multicritical potentials should be accessible from the simple cubic potential 
of (2.4). (On the other hand, it is becoming increasingly apparent that matrix models 
involving an external matrix field possess a rich multicritical structure [6, 17].) 

Finally, we complete the correspondence between the Kontsevich model and topo- 
logical gravity by noting that the familiar selection rule for correlation functions of 
scaling operators at genus g 

<IK">» = n(aJ NH. (213) 



** For related comments, see section 1.1 of ref. 2. 
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(where [• • •] means the contribution at genus g), namely 

J2n k (k-l) = 3g-3 , (2.14) 
k 

is readily extracted from results in the appendix of ref. 6. As is well known, eqn. (2.14) 
and the Lq constraint are all one needs to find the string susceptibility and the scaling 
dimensions of operators. 



3. Unitary matrices. 

The success of the analysis so far suggests that one should investigate other models 
involving an external field. One such case, considered in another context in refs. 8, 6, 
is that of a unitary matrix, 

Z = JvUexpTr(A^U + U^Aj , (3.1) 

where A is now an arbitrary matrix. As mentioned earlier, we know of no simple 
geometrical interpretation for such an integral. One approach is to write the unitary 
matrix in terms of a hermitian matrix, for example as U — e tH [18], or U = (1 + 
iH)/(l — iH) [19], but neither of these has a particularly attractive expansion in 
terms of surfaces. 

The partition function Z satisfies the differential equations 

d 2 Z 

r- = 5 ac Z . (3.2) 

dA a Mc 

After changing variables in these equations to the eigenvalues A a = fi 2 of A^A, we 
get the Schwinger- Dyson equations for Z, 

d 2 Z 1 tdZ dZ\ if dZ \ 

Guided by experience, we expect that we should factor Z along the lines of (2.6); and 
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indeed, after a small amount of trial and error, we find it useful to write 



Z = H^ a + V b r 1/2 exp(2j2vb) Y. (3.4) 

a, b b 

This factor is what one gets by expanding U about the saddle point of the action as 
the exponential of a hermitian matrix and dropping terms higher than quadratic, and 
in a sense (3.4) is the natural analog of the factorization (2.6). 

It was seen in ref. 6 that Y shares with Y the property of depending only on 
the tft's. This encourages us to try changing variables in (3.3), and we find that Y 
satisfies 



i E 2fc+2j+6 8 ^ Y + E ~^+I d k Y 
k,j=0 ^ a fc=0 /ia 

u-n V Pa Ma / 



(3.5) 



fc=0 

which yields the equations 



+ —Y = , 



n = : Y,^ k + ^ d k Y + T6 Y = -doY , 

oo n 

n>l : E( A; + i)^+- F + 3E^-i^ F = -^ F - 



(3.6) 



fc=0 fc=l 



We recognize these as Virasoro constraints, namely 

L n Y = -d n Y, n>0. (3.7) 

These equations can be put in the standard form L n Y = by a shift of to, which 
makes explicit that (3.7) corresponds to an expansion about the point to — 2, U — 
for i > 1. (Contrast this with the usual situation in hermitian matrix models, where 
to = x is a free parameter, and one other t m , for some m > 0, is fixed to a non-zero 
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value.) This is an exceedingly trivial point, the m = unitary matrix model, which 
corresponds naively to a matrix model with potential V(U) = constant (the m — 1 
model has V — U + W); nevertheless it is well defined, since the integration is over 
a compact group. 

The Virasoro constraints are a convenient starting point for extracting useful 
information about the model. For instance, we can follow the example of ref. 3 and 
derive from them recursion relations for the correlation functions of scaling operators. 
For the m = model, these take the form 

(r n IJ T k)g = ~ + i)( T j+n II Tk ^9 
kes jes k^j 

n ( ^ 1 (3.8) 

~ 2 yi \ ( T j~ iT n-j n Tfc ^-i + \ yi ^j- 1 n T k)gA r n-j n \ ■ 

j=i kes s=xuy keX keY [ 

V 9 = 91+92 ) 

(The genus dependence of these correlation functions has been inserted by hand, so as 
to correspond to the conventional 1/N expansion.) Note that the recursion relations 
are precisely sufficient to determine all the correlation functions as pure numbers - 
i.e., all the operators are redundant. In a loose sense, therefore, this model is also 
"topological." 

There is also a selection rule, 

J2kn k = g-1 <p=> J2 kt k d kHY] g = (g-l)[\ogY] g , (3.9) 
k k 

which, like (2.14), was found in ref. 6. As a check, note that (3.9) and (2.14) are 
merely different linear combinations of Lq and the dilaton equation J2k n k = % ~ 2*9 ■ 
(In general, one finds the selection rule for the m th multicritical model by eliminating 
t m between the two equations.) 

★ We thank R. Dijkgraaf for pointing this out to us. 
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An important feature of (3.7) is that because there is no L_i equation, the con- 
stant in Lq is not determined simply by the requirement of closure of the Virasoro 
algebra. In ref. 11 the constant is considered a free parameter of the solution, but 
our derivation singles out the value jq. 

There is another difference between our results and theirs, which we do not really 
understand. The constraints found here act directly on Y, whereas in ref. 11 the 
partition function is a product of two tau functions, each of which is annihilated 
separately by the L n 's. 

4. Symmetric matrices and unoriented surfaces. 

In this final section we discuss an interesting example of a model which can be 
solved only partially by the methods we have been using: the real-symmetric matrix 
model. This has been investigated in several papers [20, 21, 22, 23], and has a well- 
known geometric interpretation as the sum over unoriented surfaces. It appears that 
this model in an external field has not been looked at before, so we permit ourselves 
to stray somewhat from the subject of Virasoro constraints and discuss this case in 
a fair amount of detail. Much of the analysis runs parallel to that of the hermitian 
matrices, and wherever possible we use the same notation 

The starting point is the integral 




for some real symmetric matrix X. In this section we need to count powers of N in 
the topological expansion, so we include an N in the exponent. 

The derivation of the Schwinger-Dyson equations from eqn. (4.1) is more subtle 
than for hermitian matrices because the matrix elements of M (or X) are no longer 
independent, and when we differentiate with respect to off-diagonal elements there is 



9 



an extra factor of 2. For example, 



-Tr (XM) = \ . (4.2) 



9Mij { Xu i = j 

As a consequence, the equations of motion have a complicated form when written in 
terms of derivatives with respect to the matrix elements of X. For the same reason we 
have to be careful when we change variables to eigenvalues, which requires evaluating 
such quantities as d 2 x a /dX pq dX rs . Remarkably, the final answer is both simple and 
familiar: 

d 2 Z 1^ 1_ ( dZ_ _ dZ_\ _ 2 u 3 ) 

dx l 2^ a x a -xb\dx a dx b J g 

However, there is a major difference between this and eqns. (2.3), (3.3): namely, the 
factor of \ multiplying the sum. This small change makes an important qualitative 
difference to the genus expansion. 

Before we solve (4.3), it will prove convenient to rescale variables to A a = Agx a . 
Then, following the procedure of ref. 6, we define Z = e NF and rewrite (4.3) as an 
equation for F a = dF/d\ a * 

1 dF a ~ 2 , 1 ^ F a - F b 1 



N d\ a + Fl + 2N^ F \ a - F \ b b - 6^ Xa ■ (44) 

In the spherical limit the first term is smaller by a factor oil/N, and we can drop it. 
The resulting equation has the solution 

^-w+lM-*")' (4 - 5) 

where it is useful to introduce the notation 



b b b 

and y is found from the equation y = — 4(? 2 cr 1 ; the superscript on F a refers to the 



★ Note that F is of order N, and F a is of order 1. 
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order in 1/N . Integrating eqn. (4.5), we find that the free energy is 



y b,c 

We can attempt to solve (4.4) order by order in 1/N to find the corrections to the free 
energy, as in ref. 6. Just write F = F^ + -^F^ + • • •, insert into (4.4), and linearize. 
In the solutions of unitary and hermitian matrices [6], the 1/N terms canceled and 
the leading correction was 0(l/iV 2 ). In the present case, this cancellation does not 
occur, because of the factor of |; therefore the free energy has a contribution of order 
1/N, as one would expect from topological arguments [24]. 

The resulting linear equation for is 




The right-hand side here is more complicated than anything that arose in the unitary 
or hermitian cases. We have not been able to solve this equation, and unfortunately 
it seems probable that the solution to (4.8) cannot be written in closed form. If 
this is true, then it is impossible even to write down the equations for higher 1/N 
corrections, let alone solve them. 

On the other hand there is no obstacle to solving (4.8) perturbatively in g (actu- 
ally, as a double power series expansion in g and 1/N). Rather than using diagram- 
matic perturbation theory, we can work directly from the differential equation (4.4). 
But first, note that we are principally interested in the Kontsevich-type integral 

Z = J VMexpNTr {-\AM 2 - \gM z ) = e NF , (4.9) 
where F and F are related by F = F — Tr A 3 /12g 2 . Then one can show that 
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F a = dF/dXa satisfies 



, (4.10) 



where /i a , A a , and Z a are defined in section 2. To derive the perturbation expansion, 
it is convenient to rearrange this into the form 

F - Za 1 Ag2 ( 1 yr Fa ~ Fb i f 2 i lgFa 1 run 

a " 4/i a ^ ^2iV^A a -Ab a iV 9A 1 

With some help from Mathematica, we find that the first few terms of the solution 
at order 1/N are 



1 1 ) o 

6N 3 (i (ii/i c (Co + w)(w + + Ma) J + °' 9 ' ' 

(1,0, C / 

where we have used the notation Sk = jf Y2b "Tf- This can be written in a more concise 



form in terms of the "shifted" eigenvalues v a (= a//i 2 + y)> 



iV 4iV ^ 11 ^ 4iV2 VhVc [v h + „ c ) 

-il y + 0(/) . 

6Ar3 the VaVhV( ^ Va + ^) + + ^ 



(4.13) 



To recover (4.12) from this, use the perturbative expansions of u a , 

Va = HaU + \y^2+-'\ , (4-14) 
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and of y, 

y = -V± £ ^2^1/2 = "V ^ ~ ^ + •••)• ( 4 - 15 ) 

This latter equation can be solved iteratively to give y in terms of the s^'s. 

The notable feature of (4.13) is that it cannot be expressed in terms of just the 
traces S2fc+i> m contrast to the unitary and hermitian cases. In fact, it would appear 
that at each order in g 2 we encounter new and more complicated invariants of the 
eigenvalues. Because the solution in its concise form (4.13) is so simple, one is tempted 
to guess the general term at order g 2n . But this simplicity is misleading, and even at 
0(g 6 ) we do not know the form of the solution. The obvious ansatz, namely 

N 4 ^',l' a U b U c U d (v a + U b )---(l/ c + U d ) ' ^ 4 ' 16 ' ) 

a,b,c,d 

can be ruled out on dimensional grounds, as we need a term of order 1/z/ 9 , not 1/z/ 10 . 

We finish by asking whether it is possible to find Virasoro-type constraints for 
unoriented surfaces. We have the first ingredient, namely the Schwinger-Dyson equa- 
tions (4.4). Next we must identify how the integral (4.9) depends on the eigenvalues 
of A, and "change variables." In the previous examples we were lucky: a perturbative 
analysis showed that the integrals were functions solely of the t^s, and the change 
of variables was simple. In the present case, we don't even know the complete set of 
variables; but we do know that it is complicated, and the change of variables would 
likely present a thorny challenge. For these reasons, we think that further progress 
in this direction is impossible. 
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5. Discussion. 



In this paper we have demonstrated an efficient method of deriving Virasoro con- 
straints satisfied by the Kontsevich model, which we generalized to the unitary matrix 
model (3.1). Now we ask whether this has revealed anything about the geometrical 
interpretation of unitary matrices. To answer this, we consider first the hermitian 
case (2.4), which is well understood. 

There is an important point about the Kontsevich model which we wish to empha- 
size. The great majority of papers on matrix models involve studying them near their 
critical points, which is where the perturbation expansion diverges and the Feynman 
diagrams have an interpretation as continuum surfaces discretized by large numbers 
of infinitesimal triangles. This is not the situation in the Kontsevich model, where 
one interprets a Feynman diagram of low order in perturbation theory as a surface 
with a finite number of punctures (equal to the number of faces on the diagram). This 
is a topological picture, in which all surfaces of given genus and number of punctures 
are represented by a single "ribbon graph"; and it yields topological information, i.e., 
the intersection numbers. In addition, the integral (2.4) also possesses a conventional 
critical point, studied in ref. 6, with an interpretation in terms of continuum surfaces. 

This situation is reminiscent of the well-known Penner model, whose perturbative 
"topological" limit encodes the virtual Euler characteristic of surfaces [14, 15, 16], 
whereas the continuum limit corresponds to a c = 1 string in a compactified target 
space of a particular radius [25] . 

Now we conjecture that the unitary matrix model of section 3 is a third example 
of this type. The integral certainly possesses the two appropriate limits. From the 
"perturbative" limit we derive Virasoro constraints. The work of Hollowood et al. [11] 
gives an independent confirmation that these are the correct constraints to describe 
the multicritical unitary matrix models, which we wish to interpret as some ensemble 
of surfaces coupled to matter. The non-trivial step is to suggest that this structure 
arises here because the integral (3.1) has a representation, analogous to that of the 
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Kontsevich integral, in terms of the "moduli space" of the same (unknown) ensemble. 
However, we leave open the question of what this ensemble might be. 

It is very disappointing that we were unable to repeat the analysis for symmetric 
matrices, but the moduli space of unoriented surfaces is known to be a lot more 
complicated than that of oriented surfaces.* Nevertheless, it has been instructive to 
study this model for the subtleties it reveals in the hermitian and unitary matrix 
models. 

As a possible generalization of this work, one could try to extend these results 
to the d = —2 external field model [17]. This would not be trivial, however, as that 
model was solved by a different method. 

The original motivation for the research that led to this paper was to find a 
geometrical picture for unitary matrix models. This is a goal that has long been 
sought by many others (e.g. [12, 18, 27]). Unfortunately a satisfactory answer eludes 
us still, but we offer our insights in the hope that others can continue this work to 
completion. 

Note added: While this paper was being prepared, we learned that A. Marshakov, 
A. Mironov, and A. Morozov [28], and M. Kontsevich [29], have independently derived 
these same results for hermitian matrices. In addition, S. Kharchev, A. Marshakov, 
A. Mironov, A. Morozov, and A. Zabrodin [30] have recently proposed a generalized 
Kontsevich model to describe all multi-matrix models. 

Acknowledgements: We are grateful to the authors of ref. 1 1 for sharing their results 
with us prior to publication, and in particular to Tim Hollowood and Andrea Pasquin- 
ucci for valuable discussions. We would also like to thank Ulf Danielsson, Robbert 
Dijkgraaf, Jacques Distler, Miguel Martin-Delgado, Herman Verlinde, and Ed Witten 
for frequent help and explanations along the way. M. N. is especially indebted to Mark 
Doyle for many useful conversations, and in particular for his thorough explanation 
of the paper of Kontsevich. 



★ For example, see ref. 26. 
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